Comnient on "Gap opening in graphene by shear strain" 
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G. Cocco, E. Cadelano, and L. Colombo [Phys. Rev. B 81, 241412(R) (2010)] have suggested 
that combinations of shear and uniaxial strain can be used to open a band gap in graphene at much 
lower levels of strain than with the application of unaxial strain alone. They employed a unit cell of 
graphene in their studies and applied the Cauchy-Born rule to model external strain. Consequently, 
an important aspect of the mechanical behavior of membranes, namely buckling and wrinkling under 
external strain, and the attendant coupling with electronic structure was ignored in their analysis. 
Upon doing so, the apparent band gap that appears in the range of 15-20% shear strain under the 
Cauchy-Born assumption is shown to vanish. The gapless spectrum of graphene is found to persist 
under large shear strains as well as large combinations of shear and uniaxial strain. 

PACS numbers: 73.22.Pr, 81.05.ue, 62.25.-g 



In a recent article, Cocco et al^ have suggested that a 
combinations of shear and uniaxial strains in the range 
of 12-17% can be used to open a band gap in graphene. 
Previous studies^ indicate that uniaxial strains as high 
as 23% along the zigzag direction are required to open 
a band gap in graphene; uniaxial strains along the arm- 
chair direction do not induce a band gap. Since strains 
^ 23% are close to the failure strength of graphene,^ 
Cocco et al.^s proposal, which relies on lower levels of 
strain, presents a more practical approach for strain- 
induced band gap opening in graphene. 

In the following, I show that Cocco et a/.'s sugges- 
tion, while interesting, is unlikely to induce band gaps at 
any reasonable (i.e. well below failure) levels of strain 
in a graphene sheet. In particular, I show that nei- 
ther shear strains as high as 20% nor combinations of 
shear and uniaxial strains as high as 14% open a gap in 
graphene. These results rely on the key observation that 
shearing a graphene sheet inevitable leads to out-of-plane 
deformation of the sheet. Such bucking/ wrinkling behav- 
ior is well known from the theory of membranes^ and, 
moreover, has been observed experimentally in graphene 
membranes.^ Cocco et a/.'s studies that employ only a 
unit cell of graphene are incapable of capturing such 
wrinkling and buckling of a graphene sheet. Conse- 
quently, there are important fundamental differences be- 
tween their results and the ones I discuss next. 

The simulation cell consists of a 

49.43A X 49.43A graphene sheet (20 unit cells along 
each primitive lattice vector) consisting of 800 C atoms 
[Fig. 1]. Periodic boundary conditions are applied along 
the cell vectors. Further computational details are 
provided at the end of this Comment. An in-plane shear 
is applied to the sheet via a deformation gradient^ 
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The corresponding Lagrange strain tensor^ is 
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Two scenarios for the morphology of the deformed sheet 
are considered: (A) no relaxation of atomic positions is 
permitted, the displacement of every atom being slaved 
to the macroscopic deformation F (Cauchy-Born rule), 
(B) relaxation of atomic positions is permitted.— Sce- 
nario A is essentially identical to that of Ref.[l| apart from 
the number of unit cells in the calculation. The resulting 
morphologies for both scenarios for 7 = 0.4 are displayed 
in Fig. 1(a). Clearly, large out-of-plane displacements oc- 
cur at this level of shear. Buckled morphologies at various 
levels of shear strain are displayed in Fig. 1(b); unsurpris- 
ingly, larger applied shears lead to larger buckling ampli- 
tudes. The relative energies of the buckled morphologies 
with respect to the flat ones are displayed in Table 1. 
By buckling out-of-plane, the graphene sheet greatly re- 
duces the stored elastic energy; this is especially apparent 
at large deformations. 



TABLE I: Relative energy AE of buckled sheet with respect 
to unbuckled sheet for various applied shears 7. 
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A^; [eV/atom] 


0.1 


-0.04 


0.2 


-0.22 


0.3 


-0.67 


0.4 


-1.10 



Next, we consider the electronic structure of the flat 
and buckled sheets as a function of applied shear 7. Fig. 
2(a) displays the density of states (DOS) as a function of 
applied shear under Scenario A. As seen from this figure, 
there is no band gap opening up to 7 = 0.3. At 7 = 0.4, 
we open up a gap of ^ 1.28 eV. The DOS is qualita- 
tively similar (albeit quantitatively different) from that 




FIG. 1: (a) Morphology of the graphene sheet subjected to a shear 7 = 0.4 with (left) and without (right) relaxation of atoms 
in the simulation cell, (b) Morphologies of buckled sheets at various levels of shear employed in this study. 



of Ref.[l|. In agreement with Ref.^], we see that a band 
gap opens up in the range 7 = 0.3 — 0.4 . Fig. 2(b) dis- 
plays the DOS for the buckled sheets at the same levels of 
applied shear 7. It is immediately apparent that there is 
no longer any evidence of a band gap at 7 = 0.4. Clearly, 
there is an intimate connection between the morphology 
of the graphene sheet and its electronic structure. Con- 
straining the sheet in an artificial configuration dictated 
by the Cauchy-Born rule keeps it from attaining physi- 
cally relevant low energy states, which in turn is seen to 
lead to an artificial band gap at higher levels of shear. 

Finally for completeness, let us consider the case of 
combined shear and uniaxial strain along the armchair 
direction; application of strain along the zigzag direction 
in conjunction with shear does not induce a gap as was 
shown in Ref. 1 and verified independently in this work . 
The applied deformation gradient is chosen to be of the 
form 
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with e = 0.14 and 7 = 0.28.^^ The morphology of the re- 
laxed sheet is shown in Fig. 3(a). As seen from Fig. 3(b), 
we obtain a band gap of ~ 0.29 eV in the DOS under the 
Cauchy-Born assumption, which once again disappears 



1 


7 


0' 





1 + e 











1 



upon allowing for atomic relaxation. 

In conclusion, it appears that the results of Cocco 
et al. for gap-opening in graphene under shear, while 
strictly correct for a unit cell of graphene, are unlikely to 
be applicable in practice to a graphene sheet of reason- 
able physical dimensions. The results in this Comment 
clearly indicate that it is crucial to account for buckling 
of the graphene sheet and the intimate coupling between 
the morphology of the sheet and its electronic structure. 
Upon doing so, it would appear that opening a band gap 
in graphene once again requires extremely high strains, 
close to the failure strain of the material. Thus, the re- 
sults in this Comment, in conjunction with those drawn 
in the work of Pereira et al.^ seem to reinforce the no- 
tion that the gapless spectrum of graphene is indeed quite 
robust under mechanical deformation. 

Computational details: Electronic structure calcula- 
tions were performed with the DFTB+ code^ employing 
the Slater-Koster parameters generated by Elstner et alJ' 
A 4 X 4 X 1 Monkhorst-Pack mesh is used for Brillouin 
zone sampling. A tolerance of 10~^ Ha is used for the 
self-consistent charge cycles. All atomic positions are re- 
laxed with a force tolerance of 0.01 Ha/Bohr. The DOS 
is obtained by smearing the eigenvalue spectrum with 
Gaussians of variance a^ = 0.005. 

Useful discussions with D. Naveh are gratefully ac- 
knowledged. 
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The deformation gradient F acts on a vector X in the 
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FIG. 2: Density of states (DOS) of graphene sheets (a) with and (b) without atomic relaxation as a function of apphed shear 
7. From (a) it is clear that a band gap opens up in the range 7 = 0.3 — 0.4 if atomic relaxation is not permitted. However, as 
seen from (b) if the sheet is allowed to attain a lower energy state via out-of-plane buckling, there is no evidence of a gap at 
7 = 0.4. 
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FIG. 3: Combined shear 7 = 0.28 and uniaxial strain e = 0.14 along the armchair direction, (a) Morphology of relaxed sheet, 
(b) Density of states for unrelaxed and relaxed sheet. As seen, the band gap of ^ 0.29 eV obtained in the unrelaxed case 
vanishes upon allowing for atomic relaxation. 



reference configuration and transforms it to a vector x = 
FX in the current configuration; simultaneously, a vector 
q in the reciprocal lattice transforms as i^~ q. 
The small strain tensor of linear elasticity e is obtained 
from the Lagrange strain tensor in the limit 7^1. Since 
shear strains ~ 15-20% are certainly not small, we will have 
no further use for e. Instead, the deformation of the simula- 
tion cell and its Brillouin zone is completely characterized 
by F. For making connection with Ref.[i|, one may loosely 



interpret the shear strain (" in their notation as 7/2 here. 
Cell vectors and angles are not relaxed for now. 
It is worth noting that this strain state is close to the 
maximum that can be sustained by the sheet, as may be 
ascertained by computing the eigenvalues of the Lagrange 
strain tensor. For e = 0.15 and 7 = 0.3, which corresponds 
roughly to the case of combined uniaxial and shear strain 
of C = 15% in Ref. 1, the sheet begins to develop vacancies 
and holes. 



